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ABSTRACT 


The  phenomenon  of  string  instability  is  well  known  in  a 
platoon  of  cars  moving  forward  in  an  Automated  Highway 
System  (AHS)  .  It  is  also  known  that  ships  can  experience  a 
similar  instability  phenomenon  from  mis-coordination  of 
guidance  and  control  laws.  The  proposed  research  studies 
the  latter  phenomenon,  generalized  in  the  case  of  multiple 
marine  surface  vessels  moving  in  a  platoon.  The  guestion  of 
how  it  is  possible  for  ships  raveling  in  formation  to 
exhibit  the  phenomenon  of  string  instability  is  answered. 
Moreover,  we  examine  under  what  conditions  this  phenomenon 
can  be  exhibited,  as  well  as  how  it  can  be  prevented. 
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INTRODUCTION 


I . 


A.  MOTIVATION 

The  phenomenon  of  string  instability  is  well  known  in 
a  platoon  of  cars  moving  forward  in  an  Automated  Highway 
System  (AHS)  ,  and  it  is  also  known  that  ships  can 
experience  a  similar  instability  phenomenon.  The  main 
objective  of  an  AHS  is  to  increase  the  traffic  flow 
capacity  of  highways.  Other  benefits  are  the  increased 
safety,  the  time  savings  and  the  decrease  in  pollution.  The 
increase  of  traffic  flow  capacity  in  particular  is  an 
imperative  need  due  to  the  continuously  expanding  number  of 
vehicles  and  the  limited  capacity  of  highways  today. 

A  basic  concept  in  an  AHS  is  the  formation  of  platoons 
or  strings  of  vehicles.  A  platoon  is  a  group  of  automated 
vehicles  that  maintain  a  tight  spacing  between  themselves. 
In  a  platoon,  each  vehicle  must  safely  follow  its 
predecessor  at  a  given  distance  [1]  .  The  control  system 
must  guarantee  the  stability  of  all  vehicles  traveling 
together  in  the  platoon.  "The  stability  of  groups  of 
interconnected  systems  is  known  as  string  stability"  [1] . 
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Figure  1 .  Platoon  of  automated  vehicles  on  an 

Automated  Highway  System,  developed  by  the  PATH 
project  (From:  [2]) 

Many  researchers  have  studied  the  phenomenon  of  string 
instability  observed  in  platoons  in  the  case  of  an  AHS . 
This  phenomenon  derives  from  delays  in  communication 
between  systems,  disturbances  or  generally  mis-coordination 
phenomena.  Marine  surface  vessels  can  experience  a  similar 
instability  phenomenon  from  mis-coordination  of  guidance 
and  control  laws.  The  instability  initiates  when  either  a 
vessel  deviates  from  the  commanded  path  or  cannot  follow 
the  speed  changes  of  its  predecessor. 

Motivated  by  the  case  of  cars,  the  proposed  research 
studies  the  latter  phenomenon,  generalized  in  the  case  of 
multiple  marine  surface  vessels  moving  in  a  platoon.  The 
question  of  how  it  is  possible  for  ships  traveling  in 
formation  to  exhibit  the  phenomenon  of  string  instability 
is  answered.  Moreover,  we  examine  under  what  conditions 
this  phenomenon  can  be  exhibited  as  well  as  how  it  can  be 
prevented . 
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B .  BACKGROUND 

The  phenomenon  of  string  instability  for  multiple 
marine  surface  vessels  has  not  been  studied  before. 
Nevertheless,  it  has  been  thoroughly  studied  in  the  case  of 
AHS  [3,1,4]  .  Moreover,  there  is  previous  work  on  path 
keeping  of  separated  autonomous  underwater  vehicles 
[5,6,7].  Professor  Papoulias  has  also  analyzed  the  loss  of 
stability  for  the  case  of  one  marine  vehicle  in  the 
horizontal  plane  [8] . 

C .  OBJECTIVES 

Two  vessels  traveling  in  a  platoon  with  constant  speed 
in  a  straight  line  commanded  path  on  the  horizontal  plane 
are  considered.  One  or  both  vessels  can  deviate  from  the 
commanded  path  due  to  external  disturbances  such  as  a  big 
wave,  and  this  could  lead  to  string  instability  phenomena 
from  mis-coordination  of  guidance  and  control  laws. 

The  purpose  of  this  study  is  to  examine  these 
instabilities.  Specifically,  the  objectives  of  the  study 
are : 

a)  To  explain  a  path  control  system' s  architecture 

b)  To  introduce  a  mathematical  model  based  on  a 
proposed  autopilot  control  law  in  coordination 
with  a  guidance  law 

c)  To  carry  out  a  stability  analysis  for  the  vessels 
traveling  in  a  platoon 

d)  To  examine  under  what  conditions  the  phenomenon 
of  string  instability  can  be  exhibited  and  how  it 
can  be  prevented 
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This  thesis  contains  four  chapters.  In  Chapter  I,  the 
need  for  an  accurate  path  control  system  as  well  as  the 
architecture  of  such  a  system  is  explained.  Chapter  II 
develops  a  mathematical  formulation  of  the  problem.  The 
equations  of  motion  in  a  horizontal  plane  are  derived. 
Nomoto' s  model  is  used  to  provide  the  basis  for  an 
autopilot  control  law  that  stabilizes  the  vessel  to  a 
commanded  heading  angle  i]/Cr  capable  of  restoring  the  vessel 
to  the  commanded  path.  The  line  of  sight  guidance  scheme  is 
used  in  coordination  with  the  control  law,  and  the  combined 
guidance  and  control  law  is  derived.  Chapter  III  presents  a 
stability  analysis  for  two  different  cases:  of  one  vessel 
and  the  string  based  on  the  particular  natural  frequencies 
conr  the  damping  ratios  (  of  each  autopilot  and  the  guidance 
preview  distances  di,  d2.  Criteria  for  stability  and 
conclusions  about  the  effect  of  damping  ratio  £  and  natural 
frequency  an  are  drawn.  A  simulation  verifies  the  results. 

D.  PATH  CONTROL  SYSTEM  ARCHITECTURE 

Ship  dynamics  is  a  field  that  involves  both 
maneuverability  (controllability)  and  motion  in  waves 
(seakeeping) .  Maneuverability  is  the  capability  of  the  ship 
to  change  direction  as  a  result  of  actions  taken  by  the 
helmsman  or  the  autopilot.  Maneuvering  characteristics  of 
ships  are  complex  phenomena,  which  include  course  keeping 
and  turning  ability.  These  are  conflicting  objectives.  A 
vessel  that  is  very  stable  and  keeps  well  its  course  will 
display  a  hardness  of  turning.  On  the  other  hand,  a  vessel 
with  enhanced  turning  ability  will  be  unstable  on  course 
and  will  turn  unintentionally  due  to  environmental 
disturbances  [9]. 
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1. 


Kinds  of  Motion  Stability 


Path  course  keeping  and  motion  stability  are  closely 
related.  Generally,  we  can  classify  the  different  kinds  of 
motion  stability  of  a  marine  surface  vessel  by  examining 
the  response  on  an  external  disturbance.  In  Figure  2,  the 
vessel  is  initially  considered  to  have  a  straight-line 
motion  with  a  constant  speed  (nominal  condition)  before  a 
disturbance  is  applied.  As  it  is  implied,  the  simplest  kind 
of  motion  stability  is  the  straight-line  (dynamic) 
stability  in  which  the  vessel  continues  to  move  in  a 
straight  line  after  release  from  a  disturbance  but  at  a 
different  direction.  In  another  case,  the  directional 
stability,  the  vessel  maintains  not  only  its  straight-line 
motion  but  its  direction  too.  Finally,  in  positional  motion 
stability,  the  vessel  acquires  its  original  path  (with  the 
same  direction)  [9]. 


Figure  2.  The  kinds  of  motion  stability  (From:  [9]) 


5 


2. 


Need  for  Path  Control 


The  need  for  accurate  path  control  arises  in  naval 
applications.  Common  practice  for  ships  is  to  move  on 
prescribed  routes.  The  helmsman  normally  is  called  to 
deliver  the  ship' s  appropriate  heading  angle  and  yaw 
velocity.  Path  accuracy  is  particularly  important  when  the 
ship  has  to  operate  in  confined  waters.  Furthermore, 

small  unmanned  marine  vehicles  suitable  for  use 
in  both  naval  and  commercial  operations  have 
unique  mission  requirements  and  dynamic  response 
characteristics.  In  particular,  they  are  required 
to  be  highly  maneuverable  and  very  responsive  as 
they  operate  in  obstacle-avoidance  and  object- 
recognition  scenarios.  [10] 

Therefore,  accurate  path  control  in  confined  or  shoal 
waters,  as  well  as  in  open  sea,  is  an  imperative  need, 
which  path  control  systems  are  called  upon  to  satisfy. 

3.  System  Architecture 

A  system  with  path  control  capabilities,  in  general, 
can  be  assumed  to  consist  of  three  main  dynamic  systems 
that  are  navigation,  guidance  and  autopilot  loops,  as 
depicted  schematically  in  Figure  3.  The  navigation  loop 
processes  actual  geographical  position  and  heading 
information  (ijr)  to  provide  the  actual  path.  The  latter  is 
compared  with  the  commanded  path,  which  comes  from  path 
planning  operations,  to  create  the  positional  error.  The 
guidance  loop  is  fed  with  positional  error  and  generates  an 
appropriate  commanded  heading  angle  (i ffc)  capable  of 
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restoring  the  vessel  to  the  commanded  path.  The  autopilot 
feedback  loop  is  called  upon  to  deliver  this  commanded 
angle  [  8 ]  . 


Figure  3.  Block-diagram  representation  of  system 

architecture.  Navigation  compares  actual  and 
commanded  path.  Guidance  provides  the  appropriate 
commanded  heading  angle  ilrc.  tfr  is  the  actual  heading 
angle  (From:  [ 8 ] ) 
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II.  PROBLEM  FORMULATION 


A.  INTRODUCTION 

This  chapter  introduces  a  mathematical  model  for  the 

problem  of  path  keeping.  Equations  of  motion  for  a  marine 
surface  vessel  during  maneuvering  in  a  horizontal  plane 
will  result  in  the  basic  turning  dynamics  of  the  vessel. 
Based  on  these  dynamics,  we  suggest  a  control  law  for  the 
autopilot,  which  in  coordination  with  an  appropriate 
guidance  scheme  stabilizes  the  vessel  to  a  commanded 

heading  angle  ifrc  capable  of  restoring  the  vessel  to  the 

commanded  path.  Finally,  we  arrive  at  the  systems  of 

equations  that  describe  the  behavior  of  two  vessels  moving 
in  a  platoon  in  vector  form,  or  at  the  combined  guidance 
and  control  law. 

B.  EQUATIONS  OF  MOTION 

Maneuvering  characteristics  of  a  marine  surface  vessel 
can  be  described  by  the  equations  of  motion.  These 
equations  can  be  derived  by  Newton' s  second  law  and  are  the 
basis  for  the  autopilot  control  law. 

The  vessel  is  considered  to  be  a  rigid  body  with  only 
three  degrees  of  freedom— surge,  sway  and  yaw— and  which 
moves  in  the  horizontal  plane.  The  other  three  degrees  of 
freedom— roll ,  pitch  and  heave— are  neglected  and  not 
considered  in  this  treatment.  The  equations  can  be 
considered  either  with  reference  to  a  fixed  Cartesian 
coordinate  system  on  the  earth  or  a  fixed  system  at  the 
center  of  mass  of  the  moving  vessel,  as  illustrated  in 
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Figure  4.  Yaw  angle  ijf,  is  the  angle  between  the  ship's 
longitudinal  axis  x  and  the  fixed  axis  xo.  Furthermore,  r  is 
the  yaw  velocity  ( r  =ip  )  .  Newton's  law  expressed  in  a  fixed 
Cartesian  coordinate  system  gives 

mx  o  =  X  o 

myo  =  Yo  ( 1 ) 

Izijj  =  N  or  hr  =  N 


Xo  : 

the 

total 

force  in  xo 

direction 

Y«  : 

the 

total 

force  in  yo 

direction 

N  : 

the 

moment 

about  zq  axis 

m  : 

the 

vessel 

'  s  mass 

L  : 

the 

moment 

of  inertia 

about  zo  axis 

x»,yo  : 

acceleration  in  xo,yo 

directions 

Figure  4.  Either  a  fixed  Cartesian  coordinate  system 

on  the  earth  or  a  fixed  system  at  the  center  of 
mass  of  the  moving  vessel  can  be  used  for  deriving 
equations  of  motion 
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Equations  (1)  are  written  with  respect  to  a  fixed 
coordinate  system  at  a  particular  place.  To  dispose  of  this 
limitation,  it  is  more  convenient  to  express  them  with 
respect  to  a  fixed  system  on  the  moving  vessel.  This  will 
be  a  counterclockwise  system,  with  the  x  positive  axis  in 
the  longitudinal  direction  and  the  y  positive  axis  at  the 
starboard  of  the  vessel  (Figure  4)  .  The  previous  equations 
become  [9] 

mu  -  mvr  =  X 

mv  +  mur  =  Y  ( 2 ) 

Lr  =  N 

where  u  :  the  surge  velocity 
v :  the  sway  velocity 
r  :  the  yaw  velocity 
X  :  the  total  force  in  x  direction 
Y  :  the  total  force  in  y  direction 

Equations  (2)  can  subsequently  be  written  with  respect 
to  a  fixed  system  at  amidships  instead  of  the  center  of 
mass  of  the  vessel  [9] 

m(u  -vr  -  xGr2)  =  X 

m(v  +  ur  +  xGr)  =  Y  ( 3 ) 

Izr  +  mxG  (v  +  ur)  =  N 

where  xG  :  the  x-coordinate  of  the  center  of  mass  G. 

Total  forces  and  moments  in  Equations  (3)  ,  which  act 
during  maneuvering,  are  composed  of  fluid  forces  components 
due  to  surrounding  water,  forces  due  to  control  surfaces 
(rudders,  bow  planes,  thrusters) ,  forces  due  to  environment 
(wind,  current,  waves)  and  the  propulsion  force  ( T ) 
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(4) 


X  =  Xf+Xs+Xe+T 

y  =  yf  +  ys  +  ye 


N  =  Nf+Ns+Ne 

where  XF,YF,NF:  the  fluid  forces  components 

Xs,Ys,Ns:  the  control  surfaces  forces  components 

XEf  Ye,Ne:  the  environment  forces  components 

At  this  point,  if  we  consider  a  nominal  condition  for 
the  vessel  that  is  the  straight  line  motion  at  constant 
speed  U,  according  to  [9],  the  fluid  forces  components  can 
be  expressed  in  the  following  form 


xP  - ax- 


li  + 


dXE 


(u  -  U ) 


dll  dll 

dYF  dYF  .  dYF  dYF  . 
v  +  — —v  +  — —r  +  — —r 


dv 


dv 


dr 


dr 


(5! 


dNF  dNF  .  dNF  dNF  . 

Nf  = — -v  +  — — v  +  — —r  +  — —r 
3v  dr  dr 

and  the  control  surfaces  forces  components  are  expressed 

Xs(,5)-Xs(0)  +  di,5=0 

do 

dY 

Ys(d)  =  Ys(0)  +  -fd=Ysd 
do 

dN 

Ns(6)  =  Ns(0)  +  -fd=Ndd 
do 

where  d  :  the  rudder  angle 

U:  the  nominal  velocity  in  x  direction 


(6) 
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Assuming  no  environmental  forces  and  substituting 
Equations  (4),  (5),  (6)  into  Equation  (3),  we  arrive  at  the 
following  sway  and  yaw  linear  equations  of  motion  in  the 
horizontal  plane 

(m  -  Y-)v  -  (Y.  -  mxG)r  =  Yvv  +  (Y  -  mu)r  +  Ydd 
(Iz  -  N- )r  -  {N-  -  mxG  )v  =  Nvv  +  (Nr  -  mxGu)r  +  Ns  6 

By  nondimensionalizing  the  above  equation  in  terms  of 
water  density  p,  vessel  length  L,  and  nominal  velocity  U, 
we  take  the  dimensionless  form  that  has  the  same  identical 
form  [9] 

(m  -  Y.)v  -  (Y.  -  mxG)r  =  7v  +  (Y  -  m)r  +  YS8 
Vz-N,)r-  (N.  -  mxG )v  =  Nvv  +  (Nr  -mxG)r  +  Nd8 

C.  AUTOPILOT  CONTROL  LAW 

Sway  and  yaw  dimensionless  equations  of  motion  (8) 
yield  the  basic  turning  dynamics  of  a  marine  surface 
vessel.  Moreover,  a  control  law,  based  on  these  dynamics, 
is  needed  for  the  autopilot  so  that  vessel's  stabilization 
to  any  commanded  heading  angle  xpc  becomes  achievable. 

1 .  Turning  Dynamics 

Re-working  the  sway-yaw  Equations  (8),  the  following 
second-order  linear  differential  equations  are  obtained  [9] 

TlT2r  +  (T1+T2y  +  r  =  K6+KT3d 
Tjy  +  (Tl+T2  )v  +  V  =  KVS  +  KT48 


13 


where  TXT2  = 

T\+T2  = 


T  = 


T  = 

14 


K  = 


-K 


(Y.  -  m)(Nr  -  4 )  -  (  Yr  -  mxG  )(NV  -  mxG  ) 

Y(Nr-mxc;U)-Nv(Y-mU) 

(Y-m)(Nr-mxGU)  +  YXNr  -7  )  (Y-mxG)Nv+(Y -mU){Nv -mxG) 
Y(Nr  - mxGU) - Nv(Yr - mU )  -mxcU)-Nv(Yr - mU ) 

Yd(Nv-mxG)-N6(Yv-m) 


Y  N  -N  Y 

v  b1v 

Ys(Nr-Iz)-Ns(Yr- 

mxG ) 

Ys(Nr 

-  mxGU)  -  Ns(Yr 

-  mU ) 

ysnv-njv 

UK 

-  mxGU)  -  Nv(Yr 

-  mU ) 

Ys(Nr 

-  mxGU)  -Ns(Yr 

-  mU ) 

YV(K 

-  mxGU)  -  Nv(Yr 

-  mU ) 

The  first  of  the  above  ordinary  linear  differential 
equations  describes  apparently  the  relationship  between 
rudder  angle  5  and  yaw  velocity  r.  The  transfer  function 
between  r,  5  can  be  obtained  from  this  equation  simply  by 
writing  the  Laplace  transformation 


z(772f )  +  L{(TX  +  T2)r)  +  L(r )  =  L(Kd)  +  L{KT26) 


or 


r  K  +  KT3s 

d  ~  TxT2s2  +(Tx+T2)s  +  l 


(10) 


where  r  ,  d  :  the  Laplace  transformations  of  time-dependent 
variables . 

Performing  a  Taylor  series  expansion  in  s  and  keeping 
the  first  order  terms.  Equation  (10)  reduces  to  a 
decipherable  first  order  transfer  function  [8] 


r  b 
8  s  -a 


(11) 
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where  b  =  —  = 


K 


T  Tx+T2-T3 


ta+t2-t 


t  =  t]+t2-t3 


or  in  differential  form 

r  =  ar  +  bd  (12) 

The  last  equation  is  called  Nomoto' s  first  order  model 
and  embodies  the  radical  characteristic  of  turning  in  a 
particularly  useful  form  since  no  sway  velocity  feedback  is 
necessary . 

So  far,  we  have  concluded  that  the  vessel  turning 
dynamics  can  be  described  by  Equation  (12)  and 

xp  =  r  (13) 


2 .  Control  Law 

Autopilot  is  intended  to  deliver  the  appropriate 
commanded  heading  angle  ipcr  capable  of  restoring  the  vessel 
to  the  commanded  path.  But  the  rudder  angle  5  in  Nomoto' s 
equation  is  the  actual  rudder  angle,  which  is  not  always 
the  same  as  the  one  that  is  delivered  by  the  autopilot. 
Actually,  the  rudder  angle  5  has  a  value  5o  that  is  given 
by  a  linear  rudder  angle  control  law,  which  will  be 
discussed  later,  until  do  reaches  a  particular  limit.  Above 
this  limit,  5  is  considered  to  be  equal  to  a  saturation 
value  5sat  (Figure  5) .  This  saturation  value  is  typically  set 
at  0.4  radians . 
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5=5 


if  |  50  |  <0, 4  rad 


Of 


or 

5=5satf  if  \  50  \  >0 , 4  rad 


(14) 


5 


Figure  5.  Behavior  of  actual  rudder  angle  5.  This  is 

not  always  the  same  with  5o  that  is  delivered  by 
the  control  law 


The  above  behavior  of  rudder  angle  5  can  be 
demonstrated  simply  by  the  hyperbolic  tangent  function 

5  =  Smt  tanh(-^)  (15) 

The  linear  feedback  heading  control  law,  which  will  be 
incorporated  into  the  autopilot,  is  based  on  Equation  (12), 
and  is  of  the  form  [8] 


80  =  kt(ip  -ipc)  +  k2r 


(16) 


The 

Equations 


system' s  transfer  function 
(12),  (13),  and  (16)  as 


is 


obtained  from 


ip  -bkx 

ip c  s2  -  {a  +  bk2  )s  -  bkx 


(17) 
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This  transfer  function  stabilizes  the  vessel  to  any 


commanded  heading  angle  ipc  .  The  system  characteristic 
equation  is 

s2  ~(a  +  bk^)s  -bkl=Q  (18) 


Comparing  the  characteristic  equation  with  the 
standard  form  for  the  second-order  system 


s2  +  2  £(ons  +  (a2  =  0 

the  controller  gains  can  be  computed  as 


a  +  2tftj, 


(19) 


(20) 


The  gains  can  be  determined  from  the  desired  values  of 
con,  These  are  selected  based  on  general  guidelines  for 
second-order  system  transient  response. 


D.  EQUATIONS  FOR  GUIDANCE 

The  autopilot,  using  the  previous  control  law, 
stabilizes  the  vessel  to  an  appropriate  commanded  heading 
angle  ipc  .  This  angle  is  generated  by  a  dynamical  system 
called  guidance.  In  this  section,  we  determine  this  angle. 

The  line  of  sight  guidance  scheme  [8]  is  illustrated 
in  Figure  6.  According  to  this,  the  vessel  is  located  at 
(x,y)  at  a  distance  y  from  the  commanded  straight-line  path 
that  is  the  x-axis.  Trying  to  come  back  to  the  commanded 
path,  the  vessel  attempts  to  direct  its  longitudinal  axis 
toward  a  target  point  D  that  is  located  ahead  of  the  vessel 
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at  a  distance  d  on  the  vessel's  nominal  path.  In  this  case, 
the  commanding  heading  angle  ipc  is  the  line  of  sight  angle 


ipc  =  -  tan 


>Z 

d 


(21) 


Figure  6.  The  line  of  sight  guidance  scheme.  The 

vessel  has  deviated  from  the  commanded  path  x-axis 
and  needs  to  direct  its  longitudinal  axis  toward  a 
target  point  D  to  come  back 


Furthermore,  this  guidance  law  is  based  on  the 
inertial  deviation  rate  from  the  commanded  path  (Figure  7) 

y  =  simp  (22) 


U=1 


v| /  y  =  sin  xp 


Figure  7 .  Inertial  deviation  rate  from  the  commanded 

path 


According  to  this  guidance  law,  the  commanded  heading 
angle  is  not  constant  but  a  function  of  the  vessel's 
position  and  the  distance  d.  A  change  in  distance  d  can 
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affect  the  speed  of  guidance  law  response.  However,  the 
autopilot  has  a  limited  reaction  time  according  to  the 
specified  natural  freguency  con  and  damping  ratio  (.  Thus, 
parameter  d  must  be  chosen  properly  in  order  for  the 
desired  response  to  be  achieved. 

E.  COMBINED  GUIDANCE  AND  CONTROL  LAW 

In  this  section,  we  derive  the  combined  guidance  and 
control  law  that  describes  the  behavior  of  a  vessel  that 
has  deviated  from  the  commanded  straight-line  path  and 
tries  to  return  by  following  a  target  point  D  (Figure  6)  . 
Afterwards,  we  expand  for  two  vessels  moving  in  a  platoon. 

1.  One  Vessel 

The  control  law  (Equation  (16))  including  the  guidance 
scheme  (Equation  (21))  becomes 

c50  =  kx  (ip  +  tan-1  — )  +  k2r  (23) 

d 

The  linearized  form  of  this  equation,  for  very  small 
heading  angle  ipc  r  is 


c50  =  k^(ip  +  — )  +  k2r  (24) 

d 

Nomoto' s  Equation  (12)  combined  with  Equations  (15) 
and  (24)  gives 


r  =  ar  +  bdsat  tanh(-^-  kx (tp  +  ^)  +  k2r ) 
o .  d 


(25) 


Then,  the  system  of  equations  that  describes  the 
behavior  of  the  complete  system,  including  the  autopilot 
control  law  along  with  the  guidance  scheme,  is  given  by  [8] 
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(26) 


xp  =r 

r  =  ar  +  bdsai  tanh(-^-  kx  (yj  +  +  k2r) 

o. a 


y  =  sin  ip 

or  in  linearized  form 


ip  =  r 


r  =  ar  +  b(kl  (ip  +  — )  +  k2r) 
d 

y=M> 


It  is  desirable  to  write  the  above  system  as 
Equation  [8] 


x  =  Ax 


where  x=  [ip  ,  r  ,  y  ] 


and 


A  = 


0  1 


0 

bk .  a  +  bk2  — - 

1  2  d 


1  0 


0 


Finally,  from  Equations  (20),  it  is 

bkx  =  -con2 
a  +  bk2  =  -2£<wb 

and  matrix  A  reforms  to 


A  = 


0  1 


( o „ 


-(on  -2'Co.)n 

d 


1  0 


0 


(27) 

matrix 

(28) 

(29) 

(30) 
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2. 


Two  Vessels  in  a  Platoon 


Now  we  expand 

the  previous  concept  for 

two 

vessels . 

In 

this  case,  the 

first  vessel  attempts 

to 

direct 

its 

longitudinal  axis 

toward  a  target  point 

D, 

whereas 

the 

second  vessel  points  its  longitudinal  axis  toward  the  first 
one  (Figure  8)  .  In  this  consideration,  di  is  the  distance 
between  the  target  point  D  and  the  first  vessel,  and  the 
distance  between  the  two  vessels  always  measured  on  the 
commanded  path. 


Figure  8.  Two  vessels  in  a  platoon,  which  have 

deviated  from  the  commanded  straight  line  path  (x 
axis ) 


The  systems  of  equations  that  describe  both  vessels 

are 

Vi  =r{  Vi  =  ri 

V  =  a\r\  +  i  i  (ip  i  +^~)  +  k1^)  and  r2  =  a2r2  +  b2(kl2(ip2  +^-)  +  k22r2)  (31) 

C |  C'f  2 

ii  =V>i  y2  =^i 

or 


x  =  Ax 


(32) 
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where 


\-^P\  r  r  r  2  '  1*2  '  y 2^  r  ^  ~ 


0 


b1kn  ax  +  bxk. 
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0 

0 


0 

¥u 

dx 

0 

0 


0 

0 


0 

0 


-b2k\2 


d2 

0 


0 

0 

0 

0 


0 

0 


0 

0 


b2kX2  a2  +  b2k22 


0  0 

1  0 

b,k. 


2  12 


a2 

0 
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III.  STABILITY  ANALYSIS 


In  this  chapter,  we  analyze  the  problem  of  stability 
for  two  vessels  moving  in  a  string  with  constant  speed  on  a 
straight-line  commanded  path  (trivial  equilibrium  solution 
characterized  by  i]/=r=y=0)  .  One  or  both  vessels  can  deviate 
from  the  commanded  path  due  to  external  disturbances  such 
as  a  big  wave.  This  is  translated  to  some  initial 
conditions  others  than  the  trivial  solution.  The  system  can 
lead  to  string  instability  phenomena  from  mis-coordination 
of  guidance  and  control  laws.  Specific  criteria  for  string 
stability,  and  conclusions  about  the  effect  of  damping 
ratio  (  and  natural  frequency  con  are  drawn.  A  simulation 
verifies  the  results.  But  before  we  study  the  string,  we 
look  at  the  previous  case  of  one  vessel  trying  to  remain  on 
the  commanded  path  (Figure  6)  . 

A.  STABILITY  OF  ONE  VESSEL 


For  the  case  of  one  vessel,  we  have  already  referred 
to  the  combined  guidance  and  control  law.  Equation  (28), 
and  the  matrix  A 


A  = 


0 

1 


1 


0 


-2^„  -°^r 
d 


0 


0 


(30) 


Local  stability  properties  can  then  be  established  by  the 
eigenvalues  of  matrix  A  [8] .  The  characteristic  equation  of 
A  is 
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(33) 


A3  +  2£ft>  A2  +  co  2 A  +  —s—  =  0 


If  we  neglect  the  constant  term  — —  from  Equation 

d 

(33),  i.e.,  c/ — >  oo  or  the  guidance  law  is  eliminated. 

Equation  (33)  reduces  to  Equation  (19)  (i.e.,  the 

characteristic  equation  when  control  law  stabilizes  the 
vessel  to  any  commanded  heading  angle  ipc  without  the 
presence  of  the  guidance  law) . 

Applying  Routh' s  criterion  to  Equation  (33)  we  have 

2^co  3  -  — —  >  0  =>  d  >  — - — 

"  d  2to> 


For  d> -  eigenvalues  of  A  have  negative  real  parts  and 

2&n 

the  combined  guidance  and  control  law  provides  stability, 
i.e.,  the  vessel  follows  the  commanded  path. 

B .  STRING  STABILITY 

1 .  Analytical  Approach 

The  controller  gains  for  the  first  vessel  are  (from 
Equations  (20 ) ) 


-(flfj  +  2 ^ | ct)  ) 


(35) 


and  for  the  second  vessel 


-&> 
r  - _ 

A-i  ^ 


k  =  - 

^22 


-(@2  +  'l^>20)n^  ) 

h, 


(36) 


Using  the  previous  values  for  gains,  matrix  A  in 
Equation  (32)  turns  to  be 


A  = 


0 


1 


0 


-a)ni2  -2^con 


-co„ 


d, 


0 

0 

0 

0 


0 

0 


(0, 


0 


0 

0 

0 

0 

COn 

1 


0 

0 

0 

1 

-2^ 

0 


0 

0 

0 

0 

-co.. 


0 


(37) 


The  study  of  the  eigenvalues  of  matrix  A  can  reveal 
the  local  stability  properties  of  the  string.  Especially, 
the  characteristic  equation  of  matrix  A 

A6  +  2 .(£j%  +  £ 2colh  )A5  +  (a ^ 2  +  con2  +  4  )A4  + 

2  2 
%  (>),u 


(38) 


( — —  + — — +  2 £,coco  +2£2cu  co  )A  + 

V  d  d  51  «1  "2  "l  «27 

2  2 

(2£i%  — y—  +  2£2  co  +C0, 2)A2  + 

d2  dl 

2  2  2  2 

(0  CO  '  -  CO  CO 

/  2  ^2  Wi  2  \  i  Wi  Wt  /~v 

(co„  — —  + — —  co„  )A  +  — 1 — —  =  0 
d2  dx  ”2  dxd2 

appears  to  be  the  product  of  the  following  two  equations 


CO  2  Q)  2 

A3  +  2£,co„  A2  +  co  2A +  — ^— =  0  and  A3  +  2£,co„  A2  +  co„  2A  +  — —  =  0  (3  9) 

1  dl  2  ”2  "2  d2 
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Equations  (39)  are  apparently  of  the  form  of  Equation  (33) , 
which  is  the  characteristic  equation  for  the  case  of  one 
vessel  tryinq  to  move  in  the  commanded  path. 

The  characteristic  equation  of  A  then  becomes 

2  2 

(A3  +  2  CjO^A"  +  <wH|  2  A  +  —j— )  •  (A3  +2  £2«„2A2  +  ft>„2'A  +  — =  0  (40) 

This  conclusion  lets  us  consider  that  the  six 
eigenvalues  of  matrix  A  are  the  roots  of  Equations  (39) 
(i.e.,  the  eigenvalues  for  the  case  of  the  two  vessels 
independently  trying  to  follow  the  commanded  path) .  This 
means  that  string  stability  is  established  if  and  only  if 
both  vessels  had  stability  with  the  assumption  they  moved 
independently  (not  in  a  string,  but  separately),  i.e.: 


dx  > 


2  £>» 


d2  > 


2^ 


(41) 


2 .  Numerical  Results 


To  simplify  the  rest  of 
damping  ratio  £  to  be  identical 
conclusions  about  the  stability 
can  be  extracted  by  the  study  of 
A,  which  becomes 


the 

analysis , 

we 

select 

for 

both  autopi 

lots 

.  Safe 

properties  of 

the 

string 

the 

eigenvalues 

of 

matrix 
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A  = 


-<V  -2 


-CD,, 


d, 


0 

0 


0),, 


0 

0 


0  0  0 

0  0  0 

2 


-Q).. 


0  0 

0  0 

0  0 
1  0 


0  0 


(42) 


In  particular,  we  can  determine  the  values  of 
parameters  £  ,  t o  ,  con  ,  dx  ,  d2  for  which  eigenvalues  of  A  have 

negative  real  parts,  which  means  that  the  combined  guidance 
and  control  law  provides  asymptotic  stability.  Actually,  we 

con 

select  a  damping  ratio  £  and  a  — -  ratio  for  the 


controllers,  while  natural  frequency  t o  varies  in  a 
suitable  range.  Using  the  simulation  code  in  Appendix  A,  we 
determine,  for  each  value  of  a>H  ,  con  ,  the  minimum  ahead 

distances  dxsm ,  d2s,m  for  which  string  stability  is 

established . 


Table  1  in  Appendix  B  summarizes  the  results  for 
different  damping  ratios  and  natural  frequencies.  The  code 
has  run  using  the  following  values  for  the  parameters: 

1.  Damping  ratio  £  :  min  value  0.4,  max  1.2,  step  0.2 

2.  Natural  frequency  (o  :min  value  l,max  4.5,  step 

0.1 

Wn 

3.  Ratio  — L:  0.5,  1  or  1.5 
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Distances  d™" ,  d2,m  :  min  value  0.05,  max  2,  step 


4  . 

0.01 

Distances  d* ,  d2*  in  the  table  are  the  hypothetical  critical 
distances  if  each  vessel  moved  alone  and  not  both  in  a 
string,  and  are  to  be  compared  with  d™n ,  d2sm  .  From  Equation 
( 34 ) ,  they  are 

d*  = — - —  and  dS  = — - —  (43) 

2^ 


As  can  be  seen  from  the  results,  hypothetical  critical 
distances  d*  and  d2  are  almost  the  same  with  the  minimum 
ahead  distances  d*'m ,  d2s,m  ( the  difference  is  due  to  the  step 

of  0.01  in  variation  of  d™n  ,  d2s,m )  .  Thus,  we  conclude  again 
that  the  critical  distances  for  string  stability  are  those 
ones  as  the  vessels  move  independently  and  not  in  a  string. 
Therefore,  we  can  study  the  stability  properties  of  the 
string  by  studying  the  stability  of  each  vessel  separately. 
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3 .  Effect  of  Damping  Ratio  (  and  Natural  Frequency  an 
on  Stability 

A  quick  look  at  stability  of  the  system  is  depicted  in 


Figures 

9, 

10,  and  11, 

where  nine  families  of 

curves 

d2 

versus 

dx  were  plotted 

for 

natural  frequencies 

%  =  1 

to 

%  =  4-5 

with 

step  0.5. 

The 

first 

three  families 

are 

for 

co„ 

OJ 

ii 

o 

bi 

the 

next  three 

for 

n2  _  | 

,  and  the  last 

three 

for 

% 

(o 

n\ 

OJ 

— —  =  1 .5  . 

The 

damping  ratio  Q 

takes 

the  values  0.6, 

00 

o 

and 

% 

1.2. 
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Z0.6p0.5 


Z0.8p0.5 


30 


Zl.2p0.5 


Figure  9.  Three  stability  diagrams  for  con  =0.5  coi:  and 

£  =  0.6  ,  £  =  0.8  ,  £  =  1.2  respectively.  Stability  curves 
are  drawn  for  co  =1  to  4.5  with  step  0.5 
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Z0.6p1 


zD.8p1 
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Zl.2p1 


1.5 


D.5 


tf  =  1  .2 


=  1 


1 

1.5 

4.5 


□  5 


1.5 


dl 


Figure  10.  Three  stability  diagrams  for  con  =  coi 

£  =  0.6  ,  £  =  0.8  ,  £  =  1.2  respectively.  Stability 
are  drawn  for  a>n  =1  to  4.5  with  step  0.5 


,  and 
curves 
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z0.6p1.5 


Z0.8p1.5 
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Zl.2p1.5 


Figure  11.  Three  stability  diagrams  for  <x>n=\.5(On  and 

£  =  0.6  ,  £  =  0.8  ,  £  =  1.2  respectively.  Stability  curves 
are  drawn  for  coH|  =1  to  4.5  with  step  0.5 

By  studying  the  previous  diagrams,  we  come  to  one 
basic  conclusion:  that,  as  natural  frequency  con,  damping 
ratio  £  ,  or  both  are  increased,  the  stability  area  expands 
to  lower  values  of  di  and  d^.  We  also  observe  that  for 
values  of  natural  frequency  co  higher  than  three,  the 
stability  curves  approach  each  other,  and  thus  stability 
area  is  not  affected  significantly  for  higher  values  of  a>  . 

Figures  12  and  13  present  stability  area  of  each 
vessel  separately.  Two  families  of  curves  for  the  first  and 
the  second  vessel,  respectively  (the  second  vessel  follows 
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the  first  one  in  the  string) ,  for  different  damping  ratios 
from  0.4  to  1.2  and  step  0.2,  have  been  made.  Natural 
frequency  co  varies  in  the  range  1-4.5  and  (On  in  the  range 

0.5-2.25  (a)n  =0.5con  )  .  Figures  14  and  15  present  another 
family  of  curves  for  the  second  vessel  for  the  cases  that 
ft)  =  ft)  and  ft),,  =  1.5ft)  ,  respectively. 


Figure  12.  A  family  of  curves  for  the  first  vessel  for 

different  damping  ratios  from  0.4  to  1.2  and  step 
0.2.  Natural  frequency  a)n  varies  in  the  range 
1-4.5 
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Figure  13.  A  family  of  curves  for  the  second  vessel  for 

different  damping  ratios  from  0.4  to  1.2  and  step 
0.2.  Natural  frequency  varies  in  the  range  1- 

4.5  and  ft)  in  the  range  0.5-2.2  5  ( ft)  =0.5  ft)  )  . 

«2  — *  /7] 
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Figure  14.  A  family  of  curves  for  the  second  vessel  for 

different  damping  ratios  from  0.4  to  1.2  and  step 
0.2.  Natural  frequency  varies  in  the  range  1- 

4.5  and  to  in  the  range  1-4.5  ( a)„  =  m  )  . 

JT')  — ' '  rl'y  rlt 
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Figure  15.  A  family  of  curves  for  the  second  vessel  for 

different  damping  ratios  from  0.4  to  1.2  and  step 
0.2.  Natural  frequency  varies  in  the  range  1- 

4.5  and  at  in  the  range  1.5-6.7  5  ( eo„  =  1.5a>  )  . 

Again,  it  is  noticed  that  the  stability  area  increases 
for  higher  values  of  £  and  c o  ,  a>n  .  We  also  see  that  as 

natural  frequency  a>n  increases,  the  stability  curves 
approach  each  other,  and  thus  stability  area  is  not 
affected  significantly  for  considerably  high  values  of  a>n  . 
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4 .  Eigenvalues 


For  every  dx>  dx  and  d2>d2*r  the  system  is  stable  and 
the  vessels  will  follow  the  commanded  path.  The  combined 
guidance  and  control  law  provides  asymptotic  stability  and 
eigenvalues  of  matrix  A  (Equation  (36) )  have  negative  real 
parts.  In  the  case  that  one  at  least  of  di  and  cfc  is  smaller 
than  the  minimum  values,  the  system  becomes  unstable  and 
the  vessel  (or  vessels)  moves  with  an  oscillatory  motion  as 
a  result  of  a  complex  conjugate  pair  of  eigenvalues  with 
positive  real  parts  that  appears  in  this  case.  At  the  limit 
of  stability,  some  eigenvalues  have  zero  real  parts  and  are 
imaginary  numbers . 

As  the  distance  dx  approaches  d*  = — - —  (43)  (from 

2  £®Bl 

lower  and  higher  values) ,  one  pair  of  complex  eigenvalues 
with  positive  imaginary  parts  approaches  each  other. 

Exactly  when  dx  =  d* ,  the  eigenvalues  cross  the  imaginary 
axis  at  the  same  point.  The  eigenvalues  are  (on  i  .  Another 
pair  of  eigenvalues  with  the  opposite  imaginary  parts 
behaves  in  the  same  way.  These  eigenvalues  are  -co„i  on  the 
limit.  The  same  happens  for  the  second  vessel.  Its 
eigenvalues  will  be  ±a>„  i  . 

The  last  two  columns  of  Table  1  in  Appendix  A  show 
these  eigenvalues  at  the  limit  of  stability  when  dx  =  dx  and 
d2  =  d2  .  These  values  are  calculated  with  a  simulation  code 
and  it  is  shown  that  they  are  exactly  the  same  with  ±<x>n  i , 

±  (0  i  • 

n  2 


40 


Imaginary 

axis 


Figure  16.  Change  of  eigenvalues  as  dx^d* 

It  is  also  derived  from  the  simulation  that  the  other 
two  eigenvalues  (matrix  A  has  six  eigenvalues)  are 

A5=-2^Bi,  A6=-2^H2  (44) 

Concluding,  the  eigenvalues  of  matrix  A  are 

*1,2  =  ±Wn  i  /  *3,4  =  /  *5  =  \  =  -2£W„2  (45) 


5 .  Simulation 

The  following  are  simulation  plots  in  Simulink.  Change 
of  deviation  y,  ,  y2  from  the  commanded  path,  with  time,  has 

been  plotted  for  stability  conditions  ( dx>  d *  and  d2>  d2* )  . 
The  plots  are  for  different  initial  conditions.  Initial 
conditions  other  than  the  trivial  solution  (ij/=r=y=0)  mainly 
come  from  external  disturbances  that  cannot  be  controlled. 
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Simulation  indicates  that,  under  stability  conditions 
{di>d]*  and  d2>d2*)r  both  vessels  always  come  back  to  the 
commanded  path  after  a  brief  oscillation,  and  string 
stability  is  maintained. 


Figure  17.  Change  of  deviation  yx  ,  y2  from  the  commanded 

path  with  time  for  stability  conditions  (I.C. 
y,  =1  ,  y2  =1.5  ) 


Figure 
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Figure  19.  Change  of  deviation  yt  ,  y2  from  the  commanded 

path  with  time  for  stability  conditions  (I.C. 

y,  =  0  ,  y2  =  1 ) 


Figure  20.  Change  of  deviation  y]  ,  y2  from  the  commanded 

path  with  time  for  stability  conditions  (I.C. 
=  0.5  ,  ip2  =  0.3  ) 
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IV.  CONCLUSIONS 


This  research  considered  the  problem  of  path  keeping 
of  marine  surface  vessels  along  commanded  paths. 
Especially,  the  phenomenon  of  string  instability,  in  the 

case  of  two  vessels  moving  in  a  platoon,  was  studied. 

The  following  objectives  have  been  accomplished: 

a)  The  architecture  of  a  path  control  system  has  been 
introduced . 

b)  An  autopilot  control  law  in  coordination  with  a 

guidance  law  has  been  proposed. 

c)  A  stability  analysis  has  been  performed. 

d)  The  conditions  for  the  phenomenon  of  string 

instability  to  be  exhibited  have  been  stated. 
Moreover,  the  question  of  how  it  can  be  prevented 
is  answered. 

e)  A  simulation  was  performed  and  verified  that  under 
stability  conditions,  both  vessels  always  keep 
themselves  at  the  commanded  path  and  string 
stability  is  maintained. 

Also,  stability  analysis  has  resulted  in  the  following 
conclusions : 

a)  String  stability  is  established  if  and  only  if 

both  vessels  had  stability  when  they  moved 
independently  (not  in  a  string  but  separately) , 
and  the  critical  distances  for  string  stability 
are  those  ones  as  when  the  vessels  move 
independently,  i.e.: 
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1 


1 


(41) 


dx  > 


- ,  d2  > - 

2^„,  ^C2con 


b)  As  in  the  case  of  one  vessel,  stability  area 
expands  to  lower  values  of  di  and  ck  as  natural 
frequency  conr  damping  ratio  (,  or  both  are 

increased.  Moreover,  for  values  of  natural 
frequency  higher  than  three,  the  stability 

curves  approach  each  other,  and  thus,  stability 
area  is  not  affected  significantly  for  higher 
values  of  co„  . 


c)  The  six  eigenvalues  of  matrix  A  in  Equation  (32) 
are  the  roots  of  Equations  (39)  (i.e.,  the 

eigenvalues  for  the  case  that  the  two  vessels 
move  independently  trying  to  follow  the  commanded 
path) .  These  values  are 


2-1.2  - 


:ft>  l 


Ka  = 


--C0,,  1 


A5  =  -2£w  ,  A6  =  -2£(on_ 


(45; 


Finally,  the  main  recommendation  for  future  research 
is  as  follows: 

a)  The  stability  conclusions  of  this  work  are  with 
regards  to  an  initial  stability  or  linear 
stability.  It  remains  to  be  seen  whether  the 
primary  mechanism  of  loss  of  stability,  in  the 
form  of  bifurcations  to  periodic  solutions,  is 
affected  by  the  presence  of  multiple  vehicles  and 
their  characteristics. 
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APPENDIX  A 


%%%%%%%%%%  simulation  code  for  string  stability  area  and  plots 

oooooooooo 

9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9'9-9' 

oooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

9-9-9-9'9-9'9-9-9-9--9-9-9-9'9- 

ooooooooooooooo 

%  Christos  Angelopoulos 
%  NPS  Monterey  CA  -  June  2011 

clear 

z=0.4;  %damping  ratio 

p=l;  %wn2=p*wnl 

V= [ ] ; D2no_st= [ ] ; Dlno_st= [ ] ; Wnlno_st= [ ] ; D2st= [ ] ; Dlst= [ ] ; Wnlst= [ ] ; Wn2st= [ 
]  ; 

Wn2no_st= [ ] ; 

%calculation  of  stability  area 

for  wnl=l : 0 . 1 : 4 . 5  %natural  frequency 

wn2=p*wnl;  %natural  frequency 

for  dl=0 . 05  :  0 . 01 :  2 
for  d2=0 . 05 : 0 . 01 :  2 

A= [ 0  1000  0 ; -wnl * wnl  -2*z*wnl  - (wnl*wnl) /dl  0  0  0;1  0  0  0  0  0;0 
0  0  0  1  0;0  0  (wn2*wn2)/d2  -wn2*wn2  -2*z*wn2  -wn2 *wn2 /d2 ; 0  0  0  1  0  0]; 

%system  matrix  A 

eigA=eig (A) ; 

V= [V  max (real (eigA) ) ] ; k=length  (V) ; 
if  V(k)>0  %no  stability  case 

D2no_st= [D2no_st; d2 ] ; 

Dlno_st= [Dlno_st; dl ]  ; 

Wnlno_st= [Wnlno_st; wnl ]  ; 

Wn2no_st= [Wn2no_st;  wn2 ]  ; 
else  %  stability  case 

D2st= [D2st; d2 ] ; 

Dlst= [Dlst; dl ] ; 

Wnlst= [Wnl st; wnl ] ; 

Wn2st= [Wn2st; wn2 ] ; 
end 

end 

end 

end 

Wnl= [Wnl st  (1)  1]; 

for  i=l : length (Wnlst) -1 

if  Wnlst (i) <Wnlst (i+1) 

Wnl= [Wnl; Wnlst (i+1)  i+1];  %matrix  of  wnl  and  position  inside 
Wnlst  when  wnl  change 

if  Wnlst (i) >Wnlst (i+1) 

Wnl= [Wnl; Wnlst (i+1)  i+1]; 
else 
end 
end 

end 

Dl= [ ] ; D2= [ ] ; 

Wn2= [Wnl ( : , 1 ) *p  Wnl(:,2)]; 
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%matrix  of  dl  inside  Dlst  when  wn  change 


index=Wn2 ( : , 2 ) ; 
for  j =1 : length ( index) ; 

D1=[D1  Dlst (index (j ))]  ; 

D2=[D2  D2st (index (j ))] ; 
end 

%Results  are  Wnl , Wn2 , Dl , D2 , 

%diagrams 

dio=2 . *ones ( 1 , length (Dl ) ) ; 

dlplot= [Dl ; dio] ; 
d2plot= [ D2 ; D2 ] ; 
dlplotl= [Dl; Dl]  ; 
d2plotl= [D2 ; dio] ; 
plot(dlplot,d2plot,dlplotl,d2plotl) 
xlabel ( ' dl ’ ) 
ylabel ( ' d2 ' ) 
axis  ([0  2.2  0  2.2]) 

Z=  [  ]  ; 

for  q=l : length (Wnl ) 

Z= [ Z ; z ]  ; 

end 

dlone_vessel=l . / (2 . *Z . *Wnl ( : ,  1 ) ) ; 
d2one_vessel=l . / (2 . *Z . *Wn2 (:,1)); 
stability__matrix=  [ Z  Wnl(:,l)  Wn2(;,l)  Dl 
d2one_vessel] 

%f inal 


'  D2 '  dlone_vessel 
stability  matrix 
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APPENDIX  B 


Table  1.  Numerical  results  for:  £  0.4  to  1.2,  step  0.2.  c o 
1  to  4.5,  step  0.1.  dx  ,  d2  0.05  to  2,  step  0.01. 
(on 

— ^  :  0.5,  1  or  1.5 


£ 

% 

co„ 

n2 

d™ 

d2sim 

d* 

d2 

eigenvalues 

0.4 

1.3 

0.65 

0.97 

1.93 

0.961538 

1.923077 

± 

1.3 

± 

0.65 

0.4 

1.4 

0.7 

0.9 

1.79 

0.892857 

1.785714 

± 

1.4 

+ 

0.7 

0.4 

1.5 

0.75 

0.84 

1.67 

0.833333 

1.666667 

± 

1.5 

+ 

0.75 

0.4 

1.6 

0.8 

0.79 

1.57 

0.78125 

1.5625 

± 

1.6 

+ 

0.8 

0.4 

1.7 

0.85 

0.74 

1.48 

0.735294 

1.470588 

± 

1.7 

± 

0.85 

0.4 

1.8 

0.9 

0.7 

1.39 

0.694444 

1.388889 

± 

1.8 

± 

0.9 

0.4 

1.9 

0.95 

0.66 

1.32 

0.657895 

1.315789 

± 

1.9 

+ 

0.95 

0.4 

2 

1 

0.63 

1.26 

0.625 

1.25 

± 

2 

± 

1 

0.4 

2.1 

1.05 

0.6 

1.2 

0.595238 

1.190476 

± 

2.1 

± 

1.05 

0.4 

2.2 

1.1 

0.57 

1.14 

0.568182 

1.136364 

± 

2.2 

+ 

1.1 

0.4 

2.3 

1.15 

0.55 

1.09 

0.543478 

1.086957 

± 

2.3 

+ 

1.15 

0.4 

2.4 

1.2 

0.53 

1.05 

0.520833 

1.041667 

± 

2.4 

± 

1.2 

0.4 

2.5 

1.25 

0.5 

1.03 

0.5 

1 

± 

2.5 

± 

1.25 

0.4 

2.6 

1.3 

0.49 

0.97 

0.480769 

0.961538 

± 

2.6 

+ 

1.3 

0.4 

2.7 

1.35 

0.47 

0.93 

0.462963 

0.925926 

± 

in 

± 

1.35 

0.4 

2.8 

1.4 

0.45 

0.9 

0.446429 

0.892857 

± 

2.8 

± 

1.4 

0.4 

2.9 

1.45 

0.44 

0.87 

0.431034 

0.862069 

± 

2.9 

+ 

1.45 

0.4 

3 

1.5 

0.42 

0.84 

0.416667 

0.833333 

± 

3 

+ 

1.5 

0.4 

3.1 

1.55 

0.41 

0.81 

0.403226 

0.806452 

± 

3.1 

± 

1.55 

0.4 

3.2 

1.6 

0.4 

0.79 

0.390625 

0.78125 

± 

3.2 

+ 

1.6 

0.4 

3.3 

1.65 

0.38 

0.76 

0.378788 

0.757576 

± 

3.3 

+ 

1.65 

0.4 

3.4 

1.7 

0.37 

0.74 

0.367647 

0.735294 

± 

3.4 

± 

1.7 

0.4 

3.5 

1.75 

0.36 

0.72 

0.357143 

0.714286 

± 

3.5 

+ 

1.75 

0.4 

3.6 

1.8 

0.35 

0.7 

0.347222 

0.694444 

± 

3.6 

+ 

1.8 

0.4 

3.7 

1.85 

0.34 

0.68 

0.337838 

0.675676 

± 

3.7 

+ 

1.85 

0.4 

3.8 

1.9 

0.33 

0.66 

0.328947 

0.657895 

± 

3.8 

± 

1.9 

0.4 

3.9 

1.95 

0.33 

0.65 

0.320513 

0.641026 

± 

3.9 

+ 

1.95 

0.4 

4 

2 

0.32 

0.63 

0.3125 

0.625 

± 

4 

+ 

2 

0.4 

4.1 

2.05 

0.31 

0.61 

0.304878 

0.609756 

± 

4.1 

+ 

2.05 
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0.4 

4.2 

2.1 

0.3 

0.6 

0.297619 

0.595238 

± 

4.2 

+ 

2.1 

0.4 

4.3 

2.15 

0.3 

0.59 

0.290698 

0.581395 

± 

4.3 

± 

2.15 

0.4 

4.4 

2.2 

0.29 

0.57 

0.284091 

0.568182 

± 

4.4 

± 

2.2 

0.4 

4.5 

2.25 

0.28 

0.56 

0.277778 

0.555556 

± 

4.5 

± 

2.25 

0.4 

1 

1 

1.25 

1.26 

1.25 

1.25 

± 

1 

± 

1 

0.4 

1.1 

1.1 

1.14 

1.14 

1.136364 

1.136364 

± 

1.1 

+ 

1.1 

0.4 

1.2 

1.2 

1.05 

1.05 

1.041667 

1.041667 

± 

1.2 

+ 

1.2 

0.4 

1.3 

1.3 

0.97 

0.97 

0.961538 

0.961538 

± 

1.3 

± 

1.3 

0.4 

1.4 

1.4 

0.9 

0.9 

0.892857 

0.892857 

± 

1.4 

± 

1.4 

0.4 

1.5 

1.5 

0.84 

0.84 

0.833333 

0.833333 

± 

1.5 

+ 

1.5 

0.4 

1.6 

1.6 

0.79 

0.79 

0.78125 

0.78125 

± 

1.6 

+ 

1.6 

0.4 

1.7 

1.7 

0.74 

0.74 

0.735294 

0.735294 

± 

1.7 

± 

1.7 

0.4 

1.8 

1.8 

0.7 

0.7 

0.694444 

0.694444 

± 

1.8 

+ 

1.8 

0.4 

1.9 

1.9 

0.66 

0.66 

0.657895 

0.657895 

± 

1.9 

+ 

1.9 

0.4 

2 

2 

0.63 

0.63 

0.625 

0.625 

± 

2 

± 

2 

0.4 

2.1 

2.1 

0.6 

0.6 

0.595238 

0.595238 

± 

2.1 

± 

2.1 

0.4 

2.2 

2.2 

0.57 

0.57 

0.568182 

0.568182 

± 

2.2 

+ 

2.2 

0.4 

2.3 

2.3 

0.55 

0.55 

0.543478 

0.543478 

± 

2.3 

± 

2.3 

0.4 

2.4 

2.4 

0.53 

0.53 

0.520833 

0.520833 

± 

2.4 

± 

2.4 

0.4 

2.5 

2.5 

0.5 

0.57 

0.5 

0.5 

± 

2.5 

+ 

2.5 

0.4 

2.6 

2.6 

0.49 

0.49 

0.480769 

0.480769 

± 

2.6 

+ 

2.6 

0.4 

2.7 

2.7 

0.47 

0.47 

0.462963 

0.462963 

± 

2.7 

± 

2.7 

0.4 

2.8 

2.8 

0.45 

0.45 

0.446429 

0.446429 

± 

2.8 

± 

2.8 

0.4 

2.9 

2.9 

0.44 

0.44 

0.431034 

0.431034 

± 

2.9 

+ 

2.9 

0.4 

3 

3 

0.42 

0.42 

0.416667 

0.416667 

± 

3 

± 

3 

0.4 

3.1 

3.1 

0.41 

0.41 

0.403226 

0.403226 

± 

3.1 

± 

3.1 

0.4 

3.2 

3.2 

0.4 

0.4 

0.390625 

0.390625 

± 

3.2 

+ 

3.2 

0.4 

3.3 

3.3 

0.38 

0.38 

0.378788 

0.378788 

± 

3.3 

+ 

3.3 

0.4 

3.4 

3.4 

0.37 

0.37 

0.367647 

0.367647 

± 

3.4 

+ 

3.4 

0.4 

3.5 

3.5 

0.36 

0.36 

0.357143 

0.357143 

± 

3.5 

+ 

3.5 

0.4 

3.6 

3.6 

0.35 

0.35 

0.347222 

0.347222 

± 

3.6 

+ 

3.6 

0.4 

3.7 

3.7 

0.34 

0.34 

0.337838 

0.337838 

± 

3.7 

± 

3.7 

0.4 

3.8 

3.8 

0.33 

0.33 

0.328947 

0.328947 

± 

3.8 

± 

3.8 

0.4 

3.9 

3.9 

0.33 

0.33 

0.320513 

0.320513 

± 

3.9 

+ 

3.9 

0.4 

4 

4 

0.32 

0.32 

0.3125 

0.3125 

± 

4 

+ 

4 

0.4 

4.1 

4.1 

0.31 

0.31 

0.304878 

0.304878 

± 

4.1 

± 

4.1 

0.4 

4.2 

4.2 

0.3 

0.3 

0.297619 

0.297619 

± 

4.2 

+ 

4.2 

0.4 

4.3 

4.3 

0.3 

0.3 

0.290698 

0.290698 

± 

4.3 

+ 

4.3 

0.4 

4.4 

4.4 

0.29 

0.29 

0.284091 

0.284091 

± 

4.4 

± 

4.4 

0.4 

4.5 

4.5 

0.28 

0.28 

0.277778 

0.277778 

± 

4.5 

± 

4.5 

50 


1.136364  0.757576 

1.041667  0.694444 

0.961538  0.641026 

0.892857  0.595238 

0.833333  0.555556 

0.78125  0.520833 

0.735294  0.490196 

0,694444  0.462963 

0.657895  0.438596 

0.625  0.416667 

0.595238  0.396825 

0.568182  0.378788 

0.543478  0.362319 

0.520833  0.347222 

0.5  0.333333 

0.480769  0.320513 

0.462963  0.308642 

0.446429  0.297619 

0.431034  0.287356 

0.416667  0.277778 

0.403226  0.268817 

0.390625  0.260417 

0.378788  0.252525 

0.367647  0.245098 

0.357143  0.238095 

0.347222  0.231481 

0.337838  0.225225 

0.328947  0.219298 

0.320513  0.213675 

0.3125  0.208333 

0.304878  0.203252 

0.297619  0.198413 

0.290698  0.193798 

0.284091  0.189394 

0.277778  0.185185 


1.67  0.833333  1.666667  ± _ 1^_± _ 0J. 

1.52  0.757576  1.515152  ±  1.1  ±  0.55 


1.39  0.694444  1.388889  ± 


1.2  ± 


0.6 

1.3 

0.65 

0.65 

1.29 

0.641026 

1.282051 

± 

1.3 

± 

0.65 

0.6 

1.4 

0.7 

0.6 

1.2 

0.595238 

1.190476 

± 

1.4 

± 

0.7 

0.6 

1.5 

0.75 

0.56 

1.12 

0.555556 

1.111111 

± 

1.5 

± 

0.75 

0.6 

1.6 

0.8 

0.53 

1.05 

0.520833 

1.041667 

± 

1.6 

+ 

0.8 

0.6 

1.7 

0.85 

0.5 

0.99 

0.490196 

0.980392 

± 

1.7 

+ 

0.85 

0.6 

1.8 

0.9 

0.47 

0.93 

0.462963 

0.925926 

± 

1.8 

± 

0.9 

0.6 

1.9 

0.95 

0.44 

0.88 

0.438596 

0.877193 

± 

1.9 

± 

0.95 

0.6 

2 

1 

0.42 

0.84 

0.416667 

0.833333 

± 

2 

+ 

1 

0.6 

2.1 

1.05 

0.4 

0.8 

0.396825 

0.793651 

± 

2.1 

± 

1.05 

0.6 

2.2 

1.1 

0.38 

0.76 

0.378788 

0.757576 

± 

2.2 

± 

1.1 

0.6 

2.3 

1.15 

0.37 

0.73 

0.362319 

0.724638 

± 

2.3 

+ 

1.15 

0.6 

2.4 

1.2 

0.35 

0.7 

0.347222 

0.694444 

± 

2.4 

+ 

1.2 

0.6 

2.5 

1.25 

0.34 

0.67 

0.333333 

0.666667 

± 

2.5 

± 

1.25 

0.6 

2.6 

1.3 

0.33 

0.65 

0.320513 

0.641026 

± 

2.6 

+ 

1.3 

0.6 

2.7 

1.35 

0.31 

0.62 

0.308642 

0.617284 

± 

2.7 

+ 

1.35 

0.6 

2.8 

1.4 

0.3 

0.6 

0.297619 

0.595238 

± 

2.8 

± 

1.4 

0.6 

2.9 

1.45 

0.29 

0.58 

0.287356 

0.574713 

± 

2.9 

± 

1.45 

0.6 

3 

1.5 

0.28 

0.56 

0.277778 

0.555556 

± 

3 

+ 

1.5 

0.6 

3.1 

1.55 

0.27 

0.54 

0.268817 

0.537634 

± 

3.1 

+ 

1.55 

0.6 

3.2 

1.6 

0.27 

0.53 

0.260417 

0.520833 

± 

3.2 

± 

1.6 

0.6 

3.3 

1.65 

0.26 

0.51 

0.252525 

0.505051 

± 

3.3 

+ 

1.65 

0.6 

3.4 

1.7 

0.25 

0.5 

0.245098 

0.490196 

± 

3.4 

+ 

1.7 

0.6 

3.5 

1.75 

0.24 

0.48 

0.238095 

0.47619 

± 

3.5 

± 

1.75 

0.6 

3.6 

1.8 

0.24 

0.47 

0.231481 

0.462963 

± 

3.6 

± 

1.8 

0.6 

3.7 

1.85 

0.23 

0.46 

0.225225 

0.45045 

± 

3.7 

+ 

1.85 

0.6 

3.8 

1.9 

0.22 

0.44 

0.219298 

0.438596 

± 

3.8 

± 

1.9 

0.6 

3.9 

1.95 

0.22 

0.43 

0.213675 

0.42735 

± 

3.9 

± 

1.95 

0.6 

4 

2 

0.21 

0.42 

0.208333 

0.416667 

± 

4 

+ 

2 

0.6 

4.1 

2.05 

0.21 

0.41 

0.203252 

0.406504 

± 

4.1 

+ 

2.05 

0.6 

4.2 

2.1 

0.2 

0.4 

0.198413 

0.396825 

± 

4.2 

± 

2.1 

0.6 

4.3 

2.15 

0.2 

0.39 

0.193798 

0.387597 

± 

4.3 

± 

2.15 

0.6 

4.4 

2.2 

0.19 

0.38 

0.189394 

0.378788 

± 

4.4 

+ 

2.2 

0.6 

4.5 

2.25 

0.19 

0.38 

0.185185 

0.37037 

± 

4.5 

+ 

2.25 

0.6 

1 

1 

0.84 

0.84 

0.833333 

0.833333 

± 

1 

± 

1 

0.6 

1.1 

1.1 

0.76 

0.76 

0.757576 

0.757576 

± 

1.1 

± 

1.1 

0.6 

1.2 

1.2 

0.7 

0.7 

0.694444 

0.694444 

± 

1.2 

± 

1.2 

0.6 

1.3 

1.3 

0.65 

0.65 

0.641026 

0.641026 

± 

1.3 

+ 

1.3 

0.6 

1.4 

1.4 

0.6 

0.6 

0.595238 

0.595238 

± 

1.4 

+ 

1.4 

0.6 

1.5 

1.5 

0.56 

0.56 

0.555556 

0.555556 

± 

1.5 

± 

1.5 

0.6 

1.6 

1.6 

0.53 

0.53 

0.520833 

0.520833 

± 

1.6 

± 

1.6 

52 


0.490196  0.490196  ±  1.7  ± 


0.462963 

0.438596 


0.462963  ± 

0.438596  ± 


0.378788 

0.362319 


0.378788  ± 

0.362319  ± 


0.333333 

0.320513 


0.333333  ± 

0.320513  ± 


0.260417 

0.252525 


0.260417  ± 

0.252525  ± 


0.213675 

0.208333 


0.213675  ± 

0.208333  ± 


0.198413 

0.193798 


0.198413  ± 

0.193798  ± 


1.8  ± 
1.9  ± 


0.416667  0.416667  ± _ 2 ±_ 

0.396825  0.396825  ±  2.1  ± 


2.2  ± 
2.3  ± 


0.347222  0.347222  ±  2.4  ± 


2.5  ± 

2.6  ± 


0.308642  0.308642  ±  2.7  ± 

0.297619  0.297619  ±  2.8  ± 

0.287356  0.287356  ±  2.9  ± 

0.277778  0.277778  ± _ 3  ± 

0.268817  0.268817  ±  3.1  ± 


3.2  ± 

3.3  ± 


0.245098  0.245098  ±  3.4  ± 

0.238095  0.238095  ±  3.5  ± 

0.231481  0.231481  ±  3.6  ± 

0.225225  0.225225  ±  3.7  ± 

0.219298  0.219298  ±  3.8  ± 


3.9  ± 


0.203252  0.203252  ±  4.1  ± 


4.2  ± 

4.3  ± 


0.189394  0.189394  ±  4.4  ± 
0.185185  0.185185  ±  4.5  ± 


0.833333  0.555556  ±  1.5  ± 

0.757576  0.505051  ±  1.65  ± 

0,694444  0.462963  ±  1.8  ± 

0.641026  0.42735  ±  1.95  ± 

0.595238  0.396825  ±  2.1  ± 

0.555556  0.37037  ±  2.25  ± 

0.520833  0.347222  ±  2.4  ± 

0.490196  0.326797  ±  2.55  ± 

0.462963  0.308642  ±  2.7  ± 

0.438596  0.292398  ±  2.85  ± 

0.416667  0.277778  ±  3  ± 


0.6 

2.1 

3.15 

0.4 

0.27 

0.396825 

0.26455 

± 

3.15 

± 

2.1 

0.6 

2.2 

3.3 

0.38 

0.26 

0.378788 

0.252525 

± 

3.3 

± 

2.2 

0.6 

2.3 

3.45 

0.37 

0.25 

0.362319 

0.241546 

± 

3.45 

± 

2.3 

0.6 

2.4 

3.6 

0.35 

0.24 

0.347222 

0.231481 

± 

3.6 

± 

2.4 

0.6 

2.5 

3.75 

0.34 

0.23 

0.333333 

0.222222 

± 

3.75 

± 

2.5 

0.6 

2.6 

3.9 

0.33 

0.22 

0.320513 

0.213675 

± 

3.9 

± 

2.6 

0.6 

2.7 

4.05 

0.31 

0.21 

0.308642 

0.205761 

± 

4.05 

± 

2.7 

0.6 

2.8 

4.2 

0.3 

0.2 

0.297619 

0.198413 

± 

4.2 

+ 

2.8 

0.6 

2.9 

4.35 

0.29 

0.2 

0.287356 

0.191571 

± 

4.35 

± 

2.9 

0.6 

3 

4.5 

0.28 

0.19 

0.277778 

0.185185 

± 

4.5 

± 

3 

0.6 

3.1 

4.65 

0.27 

0.18 

0.268817 

0.179211 

± 

4.65 

± 

3.1 

0.6 

3.2 

4.8 

0.27 

0.18 

0.260417 

0.173611 

± 

4.8 

± 

3.2 

0.6 

3.3 

4.95 

0.26 

0.17 

0.252525 

0.16835 

± 

4.95 

± 

3.3 

0.6 

3.4 

5.1 

0.25 

0.17 

0.245098 

0.163399 

± 

5.1 

+ 

3.4 

0.6 

3.5 

5.25 

0.24 

0.16 

0.238095 

0.15873 

± 

5.25 

+ 

3.5 

0.6 

3.6 

5.4 

0.24 

0.16 

0.231481 

0.154321 

± 

5.4 

± 

3.6 

0.6 

3.7 

5.55 

0.23 

0.16 

0.225225 

0.15015 

± 

5.55 

± 

3.7 

0.6 

3.8 

5.7 

0.22 

0.15 

0.219298 

0.146199 

± 

5.7 

+ 

3.8 

0.6 

3.9 

5.85 

0.22 

0.15 

0.213675 

0.14245 

± 

5.85 

± 

3.9 

0.6 

4 

6 

0.21 

0.14 

0.208333 

0.138889 

± 

6 

± 

4 

0.6 

4.1 

6.15 

0.21 

0.14 

0.203252 

0.135501 

± 

6.15 

± 

4.1 

0.6 

4.2 

6.3 

0.2 

0.14 

0.198413 

0.132275 

± 

6.3 

± 

4.2 

0.6 

4.3 

6.45 

0.2 

0.13 

0.193798 

0.129199 

± 

6.45 

± 

4.3 

0.6 

4.4 

6.6 

0.19 

0.13 

0.189394 

0.126263 

± 

6.6 

± 

4.4 

0.6 

4.5 

6.75 

0.19 

0.13 

0.185185 

0.123457 

± 

6.75 

± 

4.5 

0.8 

1 

0.5 

0.63 

1.25 

0.625 

1.25 

± 

1 

+ 

0.5 

0.8 

1.1 

0.55 

0.57 

1.14 

0.568182 

1.136364 

± 

1.1 

± 

0.55 

0.8 

1.2 

0.6 

0.53 

1.05 

0.520833 

1.041667 

± 

1.2 

± 

0.6 

0.8 

1.3 

0.65 

0.49 

0.97 

0.480769 

0.961538 

± 

1.3 

± 

0.65 

0.8 

1.4 

0.7 

0.45 

0.9 

0.446429 

0.892857 

± 

1.4 

± 

0.7 

0.8 

1.5 

0.75 

0.42 

0.84 

0.416667 

0.833333 

± 

1.5 

+ 

0.75 

0.8 

1.6 

0.8 

0.4 

0.79 

0.390625 

0.78125 

± 

1.6 

+ 

0.8 

0.8 

1.7 

0.85 

0.37 

0.74 

0.367647 

0.735294 

± 

1.7 

+ 

0.85 

0.8 

1.8 

0.9 

0.35 

0.7 

0.347222 

0.694444 

± 

1.8 

+ 

0.9 

0.8 

1.9 

0.95 

0.33 

0.66 

0.328947 

0.657895 

± 

1.9 

± 

0.95 

0.8 

2 

1 

0.32 

0.63 

0.3125 

0.625 

± 

2 

+ 

1 

0.8 

2.1 

1.05 

0.3 

0.6 

0.297619 

0.595238 

± 

2.1 

+ 

1.05 

0.8 

2.2 

1.1 

0.29 

0.57 

0.284091 

0.568182 

± 

2.2 

± 

1.1 

0.8 

2.3 

1.15 

0.28 

0.55 

0.271739 

0.543478 

± 

2.3 

+ 

1.15 

0.8 

2.4 

1.2 

0.27 

0.53 

0.260417 

0.520833 

± 

2.4 

+ 

1.2 

54 


0.8 

2.5 

1.25 

0.25 

0.51 

0.25 

0.5 

± 

2.5 

+ 

1.25 

0.8 

2.6 

1.3 

0.25 

0.49 

0.240385 

0.480769 

± 

2.6 

± 

1.3 

0.8 

2.7 

1.35 

0.24 

0.47 

0.231481 

0.462963 

± 

2.7 

± 

1.35 

0.8 

2.8 

1.4 

0.23 

0.45 

0.223214 

0.446429 

± 

2.8 

+ 

1.4 

0.8 

2.9 

1.45 

0.22 

0.44 

0.215517 

0.431034 

± 

2.9 

+ 

1.45 

0.8 

3 

1.5 

0.21 

0.42 

0.208333 

0.416667 

± 

3 

± 

1.5 

0.8 

3.1 

1.55 

0.21 

0.41 

0.201613 

0.403226 

± 

3.1 

± 

1.55 

0.8 

3.2 

1.6 

0.2 

0.4 

0.195313 

0.390625 

± 

3.2 

+ 

1.6 

0.8 

3.3 

1.65 

0.19 

0.38 

0.189394 

0.378788 

± 

3.3 

± 

1.65 

0.8 

3.4 

1.7 

0.19 

0.37 

0.183824 

0.367647 

± 

3.4 

± 

1.7 

0.8 

3.5 

1.75 

0.18 

0.36 

0.178571 

0.357143 

± 

3.5 

+ 

1.75 

0.8 

3.6 

1.8 

0.18 

0.35 

0.173611 

0.347222 

± 

3.6 

+ 

1.8 

0.8 

3.7 

1.85 

0.17 

0.34 

0.168919 

0.337838 

± 

3.7 

± 

1.85 

0.8 

3.8 

1.9 

0.17 

0.33 

0.164474 

0.328947 

± 

3.8 

+ 

1.9 

0.8 

3.9 

1.95 

0.17 

0.33 

0.160256 

0.320513 

± 

3.9 

+ 

1.95 

0.8 

4 

2 

0.16 

0.32 

0.15625 

0.3125 

± 

4 

± 

2 

0.8 

4.1 

2.05 

0.16 

0.31 

0.152439 

0.304878 

± 

4.1 

± 

2.05 

0.8 

4.2 

2.1 

0.15 

0.3 

0.14881 

0.297619 

± 

4.2 

+ 

2.1 

0.8 

4.3 

2.15 

0.15 

0.3 

0.145349 

0.290698 

± 

4.3 

+ 

2.15 

0.8 

4.4 

2.2 

0.15 

0.29 

0.142045 

0.284091 

± 

4.4 

± 

2.2 

0.8 

4.5 

2.25 

0.14 

0.28 

0.138889 

0.277778 

± 

4.5 

+ 

2.25 

0.8 

1 

1 

0.63 

0.63 

0.625 

0.625 

± 

1 

+ 

1 

0.8 

1.1 

1.1 

0.57 

0.57 

0.568182 

0.568182 

± 

1.1 

+ 

1.1 

0.8 

1.2 

1.2 

0.53 

0.53 

0.520833 

0.520833 

± 

1.2 

± 

1.2 

0.8 

1.3 

1.3 

0.49 

0.49 

0.480769 

0.480769 

± 

1.3 

± 

1.3 

0.8 

1.4 

1.4 

0.45 

0.45 

0.446429 

0.446429 

± 

1.4 

+ 

1.4 

0.8 

1.5 

1.5 

0.42 

0.42 

0.416667 

0.416667 

± 

1.5 

± 

1.5 

0.8 

1.6 

1.6 

0.4 

0.4 

0.390625 

0.390625 

± 

1.6 

± 

1.6 

0.8 

1.7 

1.7 

0.37 

0.37 

0.367647 

0.367647 

± 

1.7 

+ 

1.7 

0.8 

1.8 

1.8 

0.35 

0.35 

0.347222 

0.347222 

± 

1.8 

+ 

1.8 

0.8 

1.9 

1.9 

0.33 

0.33 

0.328947 

0.328947 

± 

1.9 

± 

1.9 

0.8 

2 

2 

0.32 

0.32 

0.3125 

0.3125 

± 

2 

± 

2 

0.8 

2.1 

2.1 

0.3 

0.3 

0.297619 

0.297619 

± 

2.1 

+ 

2.1 

0.8 

2.2 

2.2 

0.29 

0.29 

0.284091 

0.284091 

± 

2.2 

± 

2.2 

0.8 

2.3 

2.3 

0.28 

0.28 

0.271739 

0.271739 

± 

2.3 

± 

2.3 

0.8 

2.4 

2.4 

0.27 

0.27 

0.260417 

0.260417 

± 

2.4 

+ 

2.4 

0.8 

2.5 

2.5 

0.25 

0.27 

0.25 

0.25 

± 

2.5 

+ 

2.5 

0.8 

2.6 

2.6 

0.25 

0.25 

0.240385 

0.240385 

± 

2.6 

± 

2.6 

0.8 

2.7 

2.7 

0.24 

0.24 

0.231481 

0.231481 

± 

2.7 

+ 

2.7 

0.8 

2.8 

2.8 

0.23 

0.23 

0.223214 

0.223214 

± 

2.8 

+ 

2.8 
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0.8 

2.9 

2.9 

0.22 

0.22 

0.215517 

0.215517 

± 

2.9 

± 

2.9 

0.8 

3 

3 

0.21 

0.21 

0.208333 

0.208333 

± 

3 

± 

3 

0.8 

3.1 

3.1 

0.21 

0.21 

0.201613 

0.201613 

± 

3.1 

± 

3.1 

0.8 

3.2 

3.2 

0.2 

0.2 

0.195313 

0.195313 

± 

3.2 

± 

3.2 

0.8 

3.3 

3.3 

0.19 

0.19 

0.189394 

0.189394 

± 

3.3 

± 

3.3 

0.8 

3.4 

3.4 

0.19 

0.19 

0.183824 

0.183824 

± 

3.4 

± 

3.4 

0.8 

3.5 

3.5 

0.18 

0.18 

0.178571 

0.178571 

± 

3.5 

± 

3.5 

0.8 

3.6 

3.6 

0.18 

0.18 

0.173611 

0.173611 

± 

3.6 

+ 

3.6 

0.8 

3.7 

3.7 

0.17 

0.17 

0.168919 

0.168919 

± 

3.7 

± 

3.7 

0.8 

3.8 

3.8 

0.17 

0.17 

0.164474 

0.164474 

± 

3.8 

± 

3.8 

0.8 

3.9 

3.9 

0.17 

0.17 

0.160256 

0.160256 

± 

3.9 

± 

3.9 

0.8 

4 

4 

0.16 

0.16 

0.15625 

0.15625 

± 

4 

± 

4 

0.8 

4.1 

4.1 

0.16 

0.16 

0.152439 

0.152439 

± 

4.1 

± 

4.1 

0.8 

4.2 

4.2 

0.15 

0.15 

0.14881 

0.14881 

± 

4.2 

+ 

4.2 

0.8 

4.3 

4.3 

0.15 

0.15 

0.145349 

0.145349 

± 

4.3 

± 

4.3 

0.8 

4.4 

4.4 

0.15 

0.15 

0.142045 

0.142045 

± 

4.4 

± 

4.4 

0.8 

4.5 

4.5 

0.14 

0.14 

0.138889 

0.138889 

± 

4.5 

± 

4.5 

0.8 

1 

1.5 

0.63 

0.42 

0.625 

0.416667 

± 

1.5 

± 

1 

0.8 

1.1 

1.65 

0.57 

0.38 

0.568182 

0.378788 

± 

1.65 

± 

1.1 

0.8 

1.2 

1.8 

0.53 

0.35 

0.520833 

0.347222 

± 

1.8 

± 

1.2 

0.8 

1.3 

1.95 

0.49 

0.33 

0.480769 

0.320513 

± 

1.95 

± 

1.3 

0.8 

1.4 

2.1 

0.45 

0.3 

0.446429 

0.297619 

± 

2.1 

± 

1.4 

0.8 

1.5 

2.25 

0.42 

0.28 

0.416667 

0.277778 

± 

2.25 

± 

1.5 

0.8 

1.6 

2.4 

0.4 

0.27 

0.390625 

0.260417 

± 

2.4 

± 

1.6 

0.8 

1.7 

2.55 

0.37 

0.25 

0.367647 

0.245098 

± 

1.7 

± 

2.55 

0.8 

1.8 

2.7 

0.35 

0.24 

0.347222 

0.231481 

± 

2.7 

± 

1.8 

0.8 

1.9 

2.85 

0.33 

0.22 

0.328947 

0.219298 

± 

2.85 

± 

1.9 

0.8 

2 

3 

0.32 

0.21 

0.3125 

0.208333 

± 

3 

± 

2 

0.8 

2.1 

3.15 

0.3 

0.2 

0.297619 

0.198413 

± 

3.15 

± 

2.1 

0.8 

2.2 

3.3 

0.29 

0.19 

0.284091 

0.189394 

± 

3.3 

± 

2.2 

0.8 

2.3 

3.45 

0.28 

0.19 

0.271739 

0.181159 

± 

3.45 

± 

2.3 

0.8 

2.4 

3.6 

0.27 

0.18 

0.260417 

0.173611 

± 

3.6 

± 

2.4 

0.8 

2.5 

3.75 

0.25 

0.17 

0.25 

0.166667 

± 

3.75 

± 

2.5 

0.8 

2.6 

3.9 

0.25 

0.17 

0.240385 

0.160256 

± 

3.9 

± 

2.6 

0.8 

2.7 

4.05 

0.24 

0.16 

0.231481 

0.154321 

± 

4.05 

± 

2.7 

0.8 

2.8 

4.2 

0.23 

0.15 

0.223214 

0.14881 

± 

4.2 

± 

2.8 

0.8 

2.9 

4.35 

0.22 

0.15 

0.215517 

0.143678 

± 

2.9 

± 

4.35 

0.8 

3 

4.5 

0.21 

0.14 

0.208333 

0.138889 

± 

3 

± 

4.5 

0.8 

3.1 

4.65 

0.21 

0.14 

0.201613 

0.134409 

± 

3.1 

± 

4.65 

0.8 

3.2 

4.8 

0.2 

0.14 

0.195313 

0.130208 

± 

3.2 

± 

4.8 
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0.8 

3.3 

4.95 

0.19 

0.13 

0.189394 

0.126263 

± 

3.3 

± 

4.95 

0.8 

3.4 

5.1 

0.19 

0.13 

0.183824 

0.122549 

± 

3.4 

± 

5.1 

0.8 

3.5 

5.25 

0.18 

0.12 

0.178571 

0.119048 

± 

3.5 

± 

5.25 

0.8 

3.6 

5.4 

0.18 

0.12 

0.173611 

0.115741 

± 

3.6 

± 

5.4 

0.8 

3.7 

5.55 

0.17 

0.12 

0.168919 

0.112613 

± 

5.55 

± 

3.7 

0.8 

3.8 

5.7 

0.17 

0.11 

0.164474 

0.109649 

± 

5.7 

± 

3.8 

0.8 

3.9 

5.85 

0.17 

0.11 

0.160256 

0.106838 

± 

5.85 

± 

3.9 

0.8 

4 

6 

0.16 

0.11 

0.15625 

0.104167 

± 

6 

± 

4 

0.8 

4.1 

6.15 

0.16 

0.11 

0.152439 

0.101626 

± 

6.15 

± 

4.1 

0.8 

4.2 

6.3 

0.15 

0.1 

0.14881 

0.099206 

± 

6.3 

± 

4.2 

0.8 

4.3 

6.45 

0.15 

0.1 

0.145349 

0.096899 

± 

6.45 

± 

4.3 

0.8 

4.4 

6.6 

0.15 

0.1 

0.142045 

0.094697 

± 

6.6 

± 

4.4 

0.8 

4.5 

6.75 

0.14 

0.1 

0.138889 

0.092593 

± 

6.75 

± 

4.5 

1 

1 

0.5 

0.5 

1.07 

0.5 

1 

± 

1 

+ 

0.5 

1 

1.1 

0.55 

0.46 

0.91 

0.454545 

0.909091 

± 

1.1 

+ 

0.55 

1 

1.2 

0.6 

0.42 

0.84 

0.416667 

0.833333 

± 

1.2 

+ 

0.6 

1 

1.3 

0.65 

0.39 

0.77 

0.384615 

0.769231 

± 

1.3 

± 

0.65 

1 

1.4 

0.7 

0.36 

0.72 

0.357143 

0.714286 

± 

1.4 

± 

0.7 

1 

1.5 

0.75 

0.34 

0.67 

0.333333 

0.666667 

± 

1.5 

± 

0.75 

1 

1.6 

0.8 

0.32 

0.63 

0.3125 

0.625 

± 

1.6 

± 

0.8 

1 

1.7 

0.85 

0.3 

0.59 

0.294118 

0.588235 

± 

1.7 

+ 

0.85 

1 

1.8 

0.9 

0.28 

0.56 

0.277778 

0.555556 

± 

1.8 

± 

0.9 

1 

1.9 

0.95 

0.27 

0.53 

0.263158 

0.526316 

± 

1.9 

+ 

0.95 

1 

2 

1 

0.25 

0.54 

0.25 

0.5 

± 

2 

± 

1 

1 

2.1 

1.05 

0.24 

0.48 

0.238095 

0.47619 

± 

2.1 

± 

1.05 

1 

2.2 

1.1 

0.23 

0.46 

0.227273 

0.454545 

± 

2.2 

± 

1.1 

1 

2.3 

1.15 

0.22 

0.44 

0.217391 

0.434783 

± 

2.3 

± 

1.15 

1 

2.4 

1.2 

0.21 

0.42 

0.208333 

0.416667 

± 

2.4 

± 

1.2 

1 

2.5 

1.25 

0.2 

0.41 

0.2 

0.4 

± 

2.5 

+ 

1.25 

1 

2.6 

1.3 

0.2 

0.39 

0.192308 

0.384615 

± 

2.6 

+ 

1.3 

1 

2.7 

1.35 

0.19 

0.38 

0.185185 

0.37037 

± 

2.7 

± 

1.35 

1 

2.8 

1.4 

0.18 

0.36 

0.178571 

0.357143 

± 

2.8 

+ 

1.4 

1 

2.9 

1.45 

0.18 

0.35 

0.172414 

0.344828 

± 

2.9 

+ 

1.45 

1 

3 

1.5 

0.17 

0.34 

0.166667 

0.333333 

± 

3 

± 

1.5 

1 

3.1 

1.55 

0.17 

0.33 

0.16129 

0.322581 

± 

3.1 

± 

1.55 

1 

3.2 

1.6 

0.16 

0.32 

0.15625 

0.3125 

± 

3.2 

+ 

1.6 

1 

3.3 

1.65 

0.16 

0.31 

0.151515 

0.30303 

± 

3.3 

+ 

1.65 

1 

3.4 

1.7 

0.15 

0.3 

0.147059 

0.294118 

± 

3.4 

± 

1.7 

1 

3.5 

1.75 

0.15 

0.29 

0.142857 

0.285714 

± 

3.5 

+ 

1.75 

1 

3.6 

1.8 

0.14 

0.28 

0.138889 

0.277778 

± 

3.6 

+ 

1.8 
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0.135135 

0.131579 

0.128205 

0.125 

0.121951 

0.119048 

0.116279 

0.113636 

0.111111 

0.5 

0.454545 

0.416667 

0.384615 

0.357143 

0.333333 

0.3125 

0.294118 

0.277778 

0.263158 

0.25 

0.238095 

0.227273 

0.217391 

0.208333 

0.2 

0.192308 

0.185185 

0.178571 

0.172414 

0.166667 

0.16129 

0.15625 

0.151515 

0.147059 

0.142857 

0.138889 

0.135135 

0.131579 

0.128205 

0.125 


0.27027  t 
0.263158  ± 

0.25641  ± 

0.25  ± 

0.243902  ± 

0.238095  ± 

0.232558  ± 

0.227273  ± 

0.222222  |  ±~ 

0.5  ± 

0.454545  ± 

0.416667  ± 

0.384615  ± 

0.357143  ± 

0.333333  ± 

0.3125  ± 

0.294118  ± 

0.277778  ± 

0.263158  ± 

0.25  ± 

0.238095  ± 

0.227273  ± 

0.217391  ± 

0.208333  ± 

0.2  ± 
0.192308  ± 

0.185185  ± 

0.178571  ± 

0.172414  ± 
0.166667  ± 

0.16129  ± 

0.15625  ± 

0.151515  ± 

0.147059  ± 

0.142857  ± 

0.138889  ± 

0.135135  ± 

0.131579  ± 

0.128205  ± 

0.125  ± 


3.7 

±  1.85 

3.8 

±  1.9 

3.9 

±  1.95 

4 

±  2 

4.1 

±  2.05 

4.2 

±  2.1 

4.3 

±  2.15 

4.4 

±  2.2 

4.5 

±  2.25 

1 

±  1 

1.1 

±  1.1 

1.2 

±  1.2 

1.3 

±  1.3 

1.4 

±  1.4 

1.5 

±  1.5 

1.6 

±  1.6 

1.7 

±  1.7 

1.8 

±  1.8 

1.9 

±  1.9 

2 

±  2 

2.1 

±  2.1 

2.2 

±  2.2 

2.3 

±  2.3 

2.4 

±  2.4 

2.5 

±  2.5 

2.6 

±  2.6 

2.7 

±  2.7 

2.8 

±  2.8 

2.9 

±  2.9 

3 

±  3 

3.1 

±  3.1 

3.2 

±  3.2 

3.3 

±  3.3 

3.4 

±  3.4 

3.5 

±  3.5 

3.6 

±  3.6 

3.7 

±  3.7 

3.8 

±  3.8 

3.9 

±  3.9 

4 

±  4 

1 

4.1 

4.1 

0.13 

0.13 

0.121951 

0.121951 

± 

4.1 

± 

4.1 

1 

4.2 

4.2 

0.12 

0.12 

0.119048 

0.119048 

± 

4.2 

± 

4.2 

1 

4.3 

4.3 

0.12 

0.12 

0.116279 

0.116279 

± 

4.3 

± 

4.3 

1 

4.4 

4.4 

0.12 

0.12 

0.113636 

0.113636 

± 

4.4 

± 

4.4 

1 

4.5 

4.5 

0.12 

0.12 

0.111111 

0.111111 

± 

4.5 

± 

4.5 

1 

1 

1.5 

0.5 

0.35 

0.5 

0.333333 

± 

1.5 

± 

1 

1 

1.1 

1.65 

0.46 

0.31 

0.454545 

0.30303 

± 

1.65 

± 

1.1 

1 

1.2 

1.8 

0.42 

0.28 

0.416667 

0.277778 

± 

1.8 

± 

1.2 

1 

1.3 

1.95 

0.39 

0.26 

0.384615 

0.25641 

± 

1.95 

± 

1.3 

1 

1.4 

2.1 

0.36 

0.24 

0.357143 

0.238095 

± 

2.1 

± 

1.4 

1 

1.5 

2.25 

0.34 

0.23 

0.333333 

0.222222 

± 

2.25 

± 

1.5 

1 

1.6 

2.4 

0.32 

0.21 

0.3125 

0.208333 

± 

2.4 

± 

1.6 

1 

1.7 

2.55 

0.3 

0.2 

0.294118 

0.196078 

± 

2.55 

± 

1.7 

1 

1.8 

2.7 

0.28 

0.19 

0.277778 

0.185185 

± 

2.7 

± 

1.8 

1 

1.9 

2.85 

0.27 

0.18 

0.263158 

0.175439 

± 

2.85 

± 

1.9 

1 

2 

3 

0.25 

0.17 

0.25 

0.166667 

± 

3 

± 

2 

1 

2.1 

3.15 

0.24 

0.16 

0.238095 

0.15873 

± 

3.15 

± 

2.1 

1 

2.2 

3.3 

0.23 

0.16 

0.227273 

0.151515 

± 

3.3 

± 

2.2 

1 

2.3 

3.45 

0.22 

0.15 

0.217391 

0.144928 

± 

3.45 

± 

2.3 

1 

2.4 

3.6 

0.21 

0.14 

0.208333 

0.138889 

± 

3.6 

± 

2.4 

1 

2.5 

3.75 

0.2 

0.14 

0.2 

0.133333 

± 

3.75 

± 

2.5 

1 

2.6 

3.9 

0.2 

0.13 

0.192308 

0.128205 

± 

3.9 

± 

2.6 

1 

2.7 

4.05 

0.19 

0.13 

0.185185 

0.123457 

± 

2.7 

± 

4.05 

1 

2.8 

4.2 

0.18 

0.12 

0.178571 

0.119048 

± 

2.8 

± 

4.2 

1 

2.9 

4.35 

0.18 

0.12 

0.172414 

0.114943 

± 

2.9 

± 

4.35 

1 

3 

4.5 

0.17 

0.12 

0.166667 

0.111111 

± 

3 

± 

4.5 

1 

3.1 
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